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ABSTRACT 

The Doppler dipole signal dominates the cosmic microwave background (CMB) anisotropy 
maps obtained by the Wilkinson Microwave Anisotropy Probe (WMAP) mission, and plays 
a key role throughout the data processing. Previously, we discovered a timing asynchronism 
of -25.6 ms between the timestamps of the spacecraft attitude and radiometer output in the 
original raw WMAP time-ordered data (TOD), which, if not corrected in following data pro- 
cessing, would generate an artificial quadrupole component (I = 2) in recovered CMB maps 
(Liu, Xiong & Li 2010). Recently, Roukema (2010b) proves that there does exist a timing- 
offset-induced error corresponding to about -25.6 ms in the WMAP calibrated TOD by study- 
ing the fluctuation variance per pixel in the temperature map recovered from the TOD as a 
function of assumed timing-offset. Here, we find evidence directly in the WMAP TOD for 
such an uncorrected timing error, possibly occurred in calculating the Doppler dipole signal 
during the WMAP team's TOD data processing. The amplitude is highly significant and is 
consistent with previous work. We also show that the uncorrected timing-offset can lead the 
WMAP CMB quadrupole to be substantially overestimated. 
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1 INTRODUCTION 

The WMAP mission makes measurements for the CMB with two 
anten nas, A and B, and records in time-order th e raw uncalibrated 
TOD ( iBennett et alJl2003d : lHinshaw et all2003l) 



draw = 9 ■ d+b. 



(1) 



With the instrument gain g and baseline b, the calibrated TOD can 
be obtained as 



d = T A ~T B +D, 



(2) 



where Ta and Tb are the antenna temperatures of the antennas A 
and B, respectively, that would be measured if the spacecraft CMB- 
frame velocity dipole were zero, D is the Doppler dipole signal 
induced by the motion of spacecraft 
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with To — 2.725 K being the CMB monopole, c the speed of light, 
v the velocity of the observer relative to the CMB rest frame, n A 
and n„ the unit direction vectors of the antenna A and B respec- 
tivelyjj The amplitude of the Doppler dipole signal is about 3 mK, 



nearly two order of magnitude stronger than the CMB anisotropy 
(~ 50 fiK). For sky map-making, the dipole signal has to be sub- 
tracted from the calibrated TOD to get the dipole-subtracted TOD 



d s = d — D . 



(4) 



The calibration parameters g and b in Eq.[TJare determ ined by 
a dipole-based calibration procedure jHinshav7et"aT]|2003h, where 
the Doppler dipole signal is initially taken as a standarqj which 
can be calculated by Eq.f3]fi'om the spacecraft direction and veloc- 
ity data. Thus an error in evaluating the dipole signal could be ex- 
pected to induce the calibration parameters in error, and then twist 
the calibrated TOD. Furthermore, a small error in evaluating the 
Doppler dipole may arouse a significant consequence on the final 
CMB temperature map via the dipole subtraction (Eq.[4ji. For exam- 
ple, in dipole calculation with Eq.fS] an antenna direction deviation 
as small as ~ 7', just about a half-pixel in the usua l WMAP res- 
olution with the resolution parameter N„ide = 512 dBennett et al.l 
l2003bh . or ~ 20 ms timing asynchronism, can cause the differen- 
tial dipole signal a nd then the dipole-subt racted TOD to be biased 
by ~ 10 - 20 fiK iUu, Xiong. & Lil2010l) . Therefore, the Doppler 
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1 The exact form of Eq. [3] should consi der the transmission imbalance, as 
given by Eq. 3 of Iffinshaw et at] 120091) . We have adopted the exact form 



in our data analysis; however, it's also confirmed that this simplified form is 
accurate enough in this work. The Sun velocity to the CMB rest frame we 
used is (-26.26, -243.71, 274.63) km/s in the Galactic coordinate. 
2 The calibration is then further improved by iteratively solve for the dipole 
and the map. 



2 H. Liu & T. P. Li 



dipole signal plays a key role throughout the WMAP data process- 
ing, and its error should be inspected very carefully. 

The WMAP spacecraft is continuously scanning the sky. For 
each datum d raw (t) or d(t) observed at a provi ded time t in the Sci - 
ence Data Table of the WMAP TOD archive dLimon et al]|2008h . 
one has to derive n A (t), n g (t) and v(t) from the Meta Data Table 
and calculate the corresponding Doppler signal by using Eq.[3] It is 
obvious that, in both data-calibration and map-making, the time t 
used to calculate the Doppler signal D(t) must be synchronous to 
the time t of the science datum d ra w(t) or d(t). But we found in 
the WMAP TOD archive that the Meta Data Table is not recorded 
simultaneously with the Science Data Table: a datum in the Meta 
Data Table is recorded 25.6ms later than the corresponding sci- 
ence datum in the Science Data Table for all bands, in other words, 
there exists a -25.6ms timing-offset between the Meta Data Ta- 
ble and the Science Data Table (Liu, Xiong & Li 2010). Such a 
timing-offset, if were incorporated into the calibration procedure 
without being corrected, would induce a significant error in the cal- 
ibrated TOD and then on the CMB map due to improperly eval- 
uating the Doppler-induced signal D(t). To check the suspected 
timing-offset-induced error in the official WMAP calibrated TOD, 
Roukema (2010b) recently studies the median per map of the tem- 
perature fluctuation variance per pixel as a function of assumed 
timing-offset and proves that the results are consistent with a - 
25.6 ms offset at lAa whereas the hypothesis of no timing error in 
the WMAP calibrated TOD is rejected at a significance of 8.5a. In 
this work, we use a more direct technique to show the existence of 
timing-offset-induced error in the official WMAP calibrated TOD. 



2 CHECKING THE TIMING OFFSET 
2.1 Method 

Let ti denotes the observation time for the z-th datum in the WMAP 
time-ordered data, d^ (ti) the released WMAP calibrated differ- 
ential datum at ti, D^ w '(ti) the real dipole component existed in 
d^ w \ti). For an assumed timing-offset At, we can compute the 
dipole component D(ti + At) at time ti + At by Eq.[3]with v, 
and n B all at ti + At and obtain the dipole-subtracted TOD 

di w) (ti;At) = d (w \ti)-D(ti+At) 

= [T A (ti) - T B (ti)] + [D {W) (U) - D(U + At)} , (5) 
and the statistic 

V(At) = J2^ W) ^ At ^ 

i 

= £){PU(ti) - T B {U)f + [D (w) (ti) - D(ti + At)] 2 + 

i 

+ 2[T A (U) - T B {ti)][D {w \ti) - D(U + At)}} . (6) 

Since the Doppler dipole signal is causally determined only by 
the motion of the spacecraft, unrelated to the astrophysical tem- 
perature fluctuations, the third term of the sum in the second 
line of Eq.[6] i.e. the cross-term between Ta(U) — T B (ti) and 
D (w) (U) - D(U + At) in Eq.|6] summed over i, is likely to vary 
randomly around zero as At is varied (we postpone further discus- 
sion of this term to Sj3]below). Without timing error in the WMAP 
data calibration, the dipole D'™' (ti) really existed in the calibrated 
TOD should be identical to that calculated with At = 0, thus we 
have D (w) (U) - D(U) = in the right side of Eq.gand V(At) is 
on average minimized at At = 0. However, if there exists an uncor- 
rected non-zero timing-offset At* between the spacecraft attitude 



Wave Band 


(At*) (ms) 


Q 


-22.33±2.52 


V 


-32.06±2.62 


w 


-21.24±2.30 


All 


-24.22±1.47 



Table 1. Timing offset in WMAP data. If there is no problem in the offi- 
cial WMAP calibrated TOD, then we should expect (At* ) = ms for all 
bands. However, this is apparently not true, indicating that there is a timing- 
asynchronism problem in the WMAP data processing. 

data and science data during the WMAP data processing, then we 
have D {w) (ti) = D(U + At*), and, consequently, V(At) is on 
average minimized at At = At* . Thus, we can produce Doppler 
signal sets with different timing-offset At and compute the statistic 
V(At) respectively. The average of all timing-offset At* that min- 
imize V(At) is a proper estimation of the amplitude of suspected 
timing error in the WMAP data processing. 

2.2 Result 

In the WMAP TOD archive, each science frame contains 15-30 
observations, and each observation last for a duration r which is 
102.4 ms, 76.8 ms, and 51.2 ms for Q-, V- and W-band, respectively. 
Following the convention used in previous works (Liu, Xiong & Li 
2010; Roukema 2010a, 2010b), the timing-offset used in this work 
is a relative one in percentage of the duration r of each observation 

At r = (t D -to)/T, (7) 

where to is the time for the instantaneous Doppler dipole signal 
D^ w \ and to the starting time of each observation. If the time used 
for the science data and spacecraft attitude data are synchronous, 
then the time for the instantaneous Doppler dipole signal D^ w > 
should be at the center of each observation and At r = 0.5 (same 
for all bands, neglecting the r difference. That's why the relative 
time At r is preferred). Similarly, if the time for the instantaneous 
Doppler dipole signal D*-™' is at the start of each observation (e.g., 
due to timing offset), then At r = 0, and in case of the end of each 
observation we have At r = 1. The relationship between At and 
At r is 

At = (At r - 0.5)r . (8) 

For red ucing the effect o f foreground emission, we use the 
KQ75 mask dGold et al.|[2009h to remove all observations with ei- 
ther side in the mask. In calculating V(At r ) with Eq.|5] Eq.[6]and 
Eq.[8] the argument At r is taken from -6 to +6 with a step of 0.1. 
For each waveband, we compute for one-day observation periods 
and record the At* that minimizes V(At r ) for each one-day pe- 
riod. For all 7-year TOD, there are about 7 x 365 At*. Finally, 
we pick out all -4 < At* < 4 (this excludes 3-7% data), 
compute the average (At*) and the standard error of (At*) by 
Y, (At* - {At*)) 2 /^N(N - 1) to give the final relative timing- 
offset estimation. The obtained (At*) are 0.282 ± 0.025, 0.083 ± 
0.034, and 0.085 ± 0.032 for Q-, V- and W-band, respectively. The 
corresponding results for the average timing-offset (At*) are listed 
in Table 1, and the histogram plots of one-year's At* is also shown 
in Fig.[T] 

In order to test the (At*) estimation, we use a secondary 
method in which we pick out all —4 < At* < 4 likewise, then 
smooth them with a 20-point-window boxcar filter, then count its 
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Figure 1. The histogram plots of one-year's averaged At* before smooth- 
ing for the WMAP Q-band (shown in the upper panel), V-band (middle 
panel), and W-band (bottom panel), respectively. 




Figure 2. Left panel: The quadrupole component of the Q-band foreground 
emission map estimated by the WMAP team, no mask is applied. Right 
panel: Same as the left panel, but the KQ75 mask has been applied. Both in 
Galactic coordinates. 

dipole-subtracted TOD d s with -25.6 ms timing-offset correction in 
evaluating the Doppler dipole D, the quadrupole decreases down to 
28.6 fiK 2 - the timing error in calibration can generate an artificial 
quadrupole signal, leading the cosmological CMB quadrupole to 
be significantly overestimated in the WMAP releasqj. 

In previous works (Liu, Xiong & Li 2010; Moss, Scott & Sig- 
urdson 2010; Roukema 2010a), it has been discovered and con- 
firmed that, given At ~ —25.6 ms and without using any CMB 
data, an artificial quadrupole component that is very similar to the 
released WMAP CMB quadrupole can be produced (Fig. [3} Q This 
strongly questions the validity of the WMAP result, and indicates 
that most of the released WMAP CMB quadrupole might be arti- 
fic ial because of the ti ming asynchronism effe c t, as demonstrated 
bv lRoukemal j2010bh ; [Liu. Xiong. & l3 j201Cl) ; lLiu & Lil ( 12011 j) 
and this work. 



histogram with 120 bins, and then fit the histogram with Gaussian 
distribution and take the center of the fitted Gaussian peak as the 
best-fit At* estimation. The obtained (At*) for Q-, V- and W-band 
from the secondary method are 0.282, 0.108 and 0.060, respec- 
tively, well consistent to the primary method. 

If there is no timing error in the official WMAP calibrated 
TOD, we should get a nearly zero (At*) for each band, but from 
Table 1 we can see that such a hypothesis is rejected at least at 8.8a 
(estimated by the band with the lowest significance for this) and 
the existence of a -25.6 ms asynchronism in the WMAP calibrated 
TOD is within 0.94<j (estimated by the average value in the last 
line of Table 1), both are consistent with lRoukemal d2010bh - 



3 CONSEQUENCE OF THE TIMING OFFSET 

The above result of checking the timing error in the WMAP TOD 
demonstrates that the official WMAP calibrated TOD d have been 
incorrectly calibrated because of incorrect usage of the Doppler 
dipole signal D. Ideally, for the uncalibrated science data d raw we 
should use t + (At*) instead of t to estimate the Doppler dipole D 
in Eq. 3 and find g and b in Eq. 1, in order to recalculate the cali- 
brated science data d rcca iib, and then derive sky maps from d reC aiib 
without any further timing offset. However, an easier, approximate 
method on large scales is to use t + (At*) instead of t in mak- 
ing maps from the (incorrectly) calibrated science data d that are 
publicly available on the WMAP website. 

To show the consequence of the timing error upon the official 
WMAP calibrated TOD, as an example, the quadrupole (power at 
I — 2) derived with pseudo-C; method from the WMAP5 offi- 
cial CMB maps is 112.7 fiK 2 (by convention, the given value is 
1(1 + l)Ci/ (2ir)), but in the modified map recovered from the new 



4 DISCUSSION 

The diagnosed timing error in this work might be due to co- 
incidental correlation between the foreground emission and the 
quadrupoles in Fig. [3] (no matter "real" or artificial). Since the 
spherical harmonics Y(l, m) with different I or m are exactly un- 
corrected, this is same to the quadrupole of the foreground emis- 
sion being correlated with Fig.(3] However, as shown by Fig. [2] the 
quadrupole component of the foreground emission is apparently 
uncorrected with any one in Fig. [3] Moreover, as shown in Table 
1, the obtained (At*) is nearly frequency-independent. These two 
facts strongly suggests that we can safely ignore the foreground 
issue in this work. 

It might also be worried that the true CMB quadrupole might 
be coincidentally correlated with the artificial quadrupole caused 
by the timing-asynchronism effect. In this case, the similarity we 
see in Fig. [3] is just something by chance. However, we have gen- 
erated 10 s randomly distributed quadrupole pairs (qi, 32) obeying 
the basic cosmic principles, especially the following two: There 
is no preferred axis, and there is no preferred spherical harmonic 
component (which means, the power expectations in fiK 2 for all 
Y(l,m) components with I — 2 and —2 ^ m ^ 2 should be 
identical). For each pair, we calculate the correlation coefficient 



3 By using zero timing-offset AT, we have obtained fully consistent 
quadrupole result to the WMAP release (using WMAP5 TOD, KQ85 mask 
and all Q, V, W bands, and such conditions are same for a none-zero AT), 
indicating that our map-making and quadrupole estimating processes are 
fine. 

4 In IMqss. S cott. & Sigurdson] i201ll) . they have obtained an artificial 
quadrupole structure that closely resembles ours, but they claim that the 
amplitude is lower than ours. 
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Figure 3. Le/? panel: The artificial CMB quadrupole component pro- 
duced by —25.6 ms timing offset. Right panel: The released WMAP CMB 
quadrupole component. Both in Galact ic coordinate and units of mK. Re- 
produced from Fig. 2 of our own work iLiu, Xiong. & Lil2010l) . 

between gi and q?, and with 10 8 pairs we see that the probabil- 
ity of coincidentally getting a correlation coefficient greater than 
0.8 (which is the correlation coefficient between the two panels in 
Fig.[3]l is about 2.5%. Thus we reject the assumption that the tim- 
ing asynchronism effect we find in this work is due to coincidental 
correlation (in other words, the cross-term in Eq.|6]( at 0.025 confi- 
dence level. 

It's also worthwhile to notice that, our work is done in the 
TOD space, which is significantly different to the temperature map 
space. Therefore, even if there were strong correlation in the tem- 
perature map space, it does not necessarily mean the same strong 
correlation in TOD space (the cross-term in Eq.[6]l. This fact could 
further decrease the worry about the correlation issues. 



5 CONCLUSION 

In previous works we have found notable systematical errors in re - 
leased WMAP temperature maps iU et all2009tlLiu & Lill2009eJ) . 
which have been confirmed by other independent works (e.g. Au- 
rich, Lustig & Steiner 2009). We then independently developed a 
self-consistent software package for WMAP data processing, and 
from the WMAP TOD produced new CMB maps which are signif- 
icantly different from the official maps (Liu & Li 2009b, 2010a). 
Our pipeline codes are already publicly released on the website of 
Tsinghua Center for Astrophysics and the CosmoCoffee forurr[__ 
Later, in searching for the source of the difference between our and 
official maps, we discovered a -25.6ms asynchronism between the 
spacecraft attitude and radiometer output timestamps in the official 
WMAP Meta and Science Data Tables, respectively; such a timing- 
offset, if not be properly corrected in data processing, should gen- 
erate serious consequences in the recovered CMB map and power 
spectrum (Liu, Xiong & Li 2010). We artificially introduced the 
-25.6 ms asynchronism into our pipeline to simulate the WMAP 
manner and then obtained fully consistent result to the WMAP 
team indeed. 

According to the well consistency between the -25.6 ms 
timing-offset directly ob served from TOD dLiu. Xiong. & Lil2010l) 
and indirectly probed bv lRoukemal fcOlObh and this work, the most 
natural explanation should be existence of an unwanted timing er- 
ror in the WMAP data. However, such an unwanted timing error 
is a lso expect e d to in troduce a blurring effect in recovered maps, 
but lRoukemal feOlOal) did not find such effect in WMAP official 

5 Our software used for WMAP data processing has been thoroughly tested 
by us, and the result has been confirmed by other works as well (Moss, 
Scott & Sigurdson 2010; Roukema 2010a). The software can be found at 
http://dpc.aire.org.cn/data/wmap/09072731/release_vl/source_code/vl/ or 
http://cosmocoffee.info/viewtopic.php?p=4525#4525. 



maps. Recently, by comparing the median per map of the fluc- 
tuation variance per pix in the tempera ture map for different as- 
sumed timing-offsets. lRoukemal l l2010bh proved that there does ex- 
ist an about -25.6 ms timing-offset in the WMAP calibrated TOD, 
which is confirmed by us in this work directly in the WMAP TOD, 
and we further show that the uncorrected timing error occurred 
at least in calculating the Doppler dipole signal. A natural expla- 
nation for above findings is that the timing-offset-induced error 
in direction did not have an effect in the compilation of the cal- 
ibrated TO D into sky maps by the actual Jupiter pointing mea- 
surements jffinshaw et"aT]|2003l) . but the calibration error, which 
on large scales consists of a timing-offset-induced pseudo-dipole 
signal already present in the calibrated TOD, remains present in 
the sky maps. 

It has to be noticed that th e error in sky-map-based determina- 
tions of CMB dipole direction dBennett et al.l2003b1 ; lrlinshaw et al.l 
120091) can also contribute to the diagnosed timing-offset. In other 
words, the timing-offset detected in this work is a synthesis of er- 
rors in timing and in dipole direction. Besides the detected timing 
offset, other observational reasons, i.e. the sidelobe pickup contam- 
ination from dipole, can also generate artificial quadrupole aligned 
with what observed in th e official CMB map , needed to be further 
removed by model fitting .Liu & Lil201 1 alibi). After template-based 
removal of artificial quadrupole, the remaining quadrupole power 
can be as low as 10.4 fiK 2 , significantly lower than 28.6 /iK 2 de- 
rived in this work, indicating that the sidelobe-pickup-induced arti- 
fact can not be ignored. The timing asynchronism in the WMAP 
raw data might be a problem special for the WMAP mission, 
whereas errors in dipole direction and sidelobe contamination are 
common problems for all CMB missions. To quantitatively esti- 
mate the effect of producing artificial CMB large-scale anisotropies 
from each possible source by thoroughly rechecking the WMAP 
data processing process is very important, not only for the WMAP- 
based cosmological study, but also for the Planck and other future 
CMB missions. 
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